Colloidal Lattice Shearing and Rupturing with a Driven Line of Particles 
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We examine the dynamics of two-dimensional colloidal systems using numerical simulations of a 
system with a drive applied to a thin region in the middle of the sample to produce a local shear. For 
a monodisperse colloidal assembly, we find a well defined decoupling transition separating a regime 
of elastic motion from a plastic phase where the particles in the driven region break away or decouple 
from the particles in the bulk, producing a shear band. For a bidisperse assembly, we find that the 
onset of a bulk disordering transition coincides with the broadening of the shear band. We identify 
several distinct dynamical regimes that are correlated with features in the velocity-force curves. 
As a function of bidispersity, the decoupling force shows a nonmonotonic behavior associated with 
features in the noise fluctuations, power spectra, and bulk velocity profiles. When pinning is added 
in the bulk, we find that the shear band regions can become more localized, causing a decoupling of 
the driven particles from the bulk particles. For a system with thermal noise and no pinning, the 
shear band region becomes more extended and the average velocity of the driven particles drops at 
the thermal disordering transition of the bulk system. 

PACS numbers: 82.70.Dd,83.80.Hj 



I. INTRODUCTION 

As the ability to manipulate small particles or groups 
of small particles has advanced, there have been a grow- 
ing number of studies examining the effect of local pertur- 
bations such as driving a single probe particle through a 
particle assembly in coUoidaii"— and granular—"— me- 
dia. The fluctuations and transport characteristics of 
the probe particle can undergo pronounced changes de- 
pending on the ordering of the surrounding media, the 
particle-particle interactions, the temperature, and the 
magnitude of the external force on the probe particle. 

In systems with short range interactions such as disks 
or granular particlesi^"— , the velocity of the probe par- 
ticle can drop under even a small change in the bulk me- 
dia density as the jamming transition is approached. For 
bulk densities below jamming, the external force needed 
to move the probe particle is small or absent and the 
probe particle velocity distribution is bimodali^. Close 
to jamming, the external force needed to push the probe 
particle through the medium becomes finite and rapidly 
increases, while the velocity fluctuations are intermittent 
and have a power law distribution that has been inter- 
preted as arising due to criticality associated with the 
jamming transition^ii^iiSii^iiS. 

For systems with longer range interactions such as 
charged colloids, there is an elastic or coupled flow regime 
at low drives where the probe particle drags the entire as- 
sembly of particles without any local tearing^. At higher 
drives there can be a well defined transition to a plastic 
response regime where the probe particle induces tearing 
rearrangements in the surrounding media, and the veloc- 
ity of the probe particle differs from that of the particles 
it is dragging. The velocity-force curves in these systems 
can be used to identify the critical driving force separat- 



ing the elastic and plastic regimes. The velocity increases 
as a power law with increasing external force in the plas- 
tic regime, similar to the behavior observed for plastic 
depinning of vortices^S and colloids^ in the presence of 
quenched disorder. 

It is also possible to consider the motion of a forced 
particle through a surrounding medium in the pres- 
ence of quenched disorder. Such a situation has 
been studied for dragging individual vortices in type- 
II superconductors2ir— , where a driven vortex interacts 
with the intrinsic pinning in the sample as well as with 
the surrounding vortices. In this case, the probe particle 
can be indirectly pinned due to its interaction with bulk 
particles that are strongly directly pinned by the sub- 
strate, causing the appearance of a threshold depinning 
force for the probe particle^. At higher drives, the bulk 
pinning can counterintuitively lower the effective drag 
on the probe particle. This occurs since the drag origi- 
nates when having a portion of the surrounding particles 
move with the probe particle, so that when the quenched 
disorder strongly immobilizes the surrounding particles, 
they can no longer move with the probe particle, leading 
to an effective decoupling of the probe particle from the 
surrounding mediaiS.. For probe particles driven through 
crystalline structures in the absence of quenched disor- 
der, shear thinning^ and directional locking effects occur 
where the amount of drag on the probe particle depends 
on the orientation of the drive with respect to the sym- 
metry directions of the crystalline structure. Flow along 
a symmetry direction produces a lower drag since the 
probe particle can move easily between particles in the 
background lattice without constantly generating new 
topological defects, while for other driving directions, the 
probe particle induces the formation of local topological 
defects, leading to increased drag2^. It was also shown 
that for driving along symmetry directions, the drag in- 
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creases as the bulk melting temperature is approached 
due to the formation of topological defects in the region 
surrounding the probe particle, leading to a local melting 
transition^. 

Here, instead of a single driven probe particle, we con- 
sider the case of a localized quasi-one-dimensional region 
of particles driven externally along a line. The drive is 
parallel to the orientation of the line. Such a geome- 
try could be created using a laser beam directed along 
the sample edge, with magnetic particles driven by a 
magnetic strip, or with particles coupled to a mechan- 
ical external drive or driven with one-dimensional ar- 
rays of optical traps^^"— . Driven line geometries have 
already been realized for particles with Yukawa interac- 
tions in dusty plasmas, where a laser beam pushes par- 
ticles only along a line^Sr— . It is also possible to create 
systems with bulk pinning and an easy flow channel in 
confined geometries^"—. A key difference between the 
point probe particle and the line probe considered here 
is that the line probe more directly mimics a shear re- 
sponse. If the system is strongly coupled, then an elastic 
response will occur in which the particles in the system 
move along with the particles in the driven line. Here we 
study the case of a line of particles driven along a sym- 
metry direction of the background lattice for monodis- 
perse ordered, bidisperse ordered, and bidisperse disor- 
dered systems, both for pin-free samples and for samples 
that have pinning present in the bulk. 
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FIG. 1: Illustration of the sample geometry for the monodis- 
perse case. Black dots: Colloids interacting via a pairwise 
Yukawa potential to form a triangular lattice. Arrows indi- 
cate the region within which an external drive in the positive 
X direction is applied to the particles. There are periodic 
boundary conditions in the x and y-directions. 



II. SIMULATION AND SYSTEM 

In Fig. [T] we illustrate the geometry of the sample, 
which consists of a two-dimensional assembly of Nc col- 
loids with periodic boundary conditions in the x and y- 
directions. The arrows in the middle of Fig. [T] indicate 
the region within which the colloids experience an ex- 
ternal drive = F^x. Colloids outside this region 
are undriven and have = 0. The colloid density is 
p = Nc/L?', where L is the size of the simulation cell. 
The particles interact via a Yukawa or screened Coulomb 
potential. In some cases, we also consider the effect of 
pinning sites that we model as localized parabolic traps 
of radius Vp. 

The dynamics of particle i arise from integrating the 
following equation of motion: 

Here rj is the damping constant, Ri(j) is the position of 
particle and Rij = |Ri — Rj|. The particle-particle 
interaction potential is V{Rij) — q^EQexp{—KRij)/Rij. 
Here Eq = Z*"^ /Aneeoao, where q is the dimensionless 
interaction strength, Z* is the effective charge of the 
colloid, e is the solvent dielectric constant, and 1/k is 
the screening length. Lengths are measured in units of 
ao, time in units of r = r]/EQ, and forces in units of 



Fq = Eo/ao. The driving force F^ — F<jx is applied 
only to particles in a region of width d at the center 
of the sample. We increase the external drive in incre- 
ments of SF]j, then hold the drive at a constant value 
for a fixed period of time and measure the average ve- 
locity of the particles within the driven channel Vc and 
of the particles in the bulk Vb- The pinning force arises 
from Np non-overlapping parabolic pinning sites that are 
placed outside the driven channel. The pinning force 
has the form Ff = Fp{Rik / Rp)Q{Rp — Rik)'R-ik, where 
Rik = |Ri — Rfel is the distance between particle i and 
the center of pinning site k, and R^fc = (R^ — 'Rk)/Rik- 
Here Rp is the pinning site radius, Fp is the maximum 
force of the pinning site, and O is the Heaviside step func- 
tion. The effects of thermal fluctuations come from the 
Langevin noise term F'^ with the properties {F'^{t)) = 
and {Ff{t)Fj'{t')) = 2-qkB5rj5{t - t'), where ks is the 
Boltzmann constant. The initial particle configurations 
are obtained by placing the colloids in a triangular ar- 
rangement, and unless otherwise noted, the average lat- 
tice constant is a = 2.0. 



III. MONODISPERSE SYSTEM 

In Fig. \^Si) we plot the average velocity of the particles 
in the driven channel region (14) versus external drive 
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FIG. 2: A monodisperse system of colloids with varied in- 
teraction strength q and with a = 2. (a) The velocities of the 
particles in the driven channel region (Vc) vs Fd, where the 
velocities are normalized by the number of particles in the 
driven region, (b) The velocities of the particles in the bulk 
outside of the driven region (14) vs Fd, where the velocities 
are normalized by the number of particles in the non-driven 
region. The curves are for q values (from left to right) q — 0.1, 
0.3, 0.5, 0.7, 0.8, 1.0, 1.2, 1.4, and 1.6. There are two regimes. 
The first is an elastic region with Vc = Vt where all the par- 
ticles move together. This is followed at higher Fd by a de- 
coupled regime where the particles in the driven line move 
past the bulk particles, while the bulk particles remain locked 
with each other. The inset in (a) shows that the driving force 
Fc at which the decoupling transition occurs increases with 
increasing q. 

Fd, and in Fig. [H^b) we plot the corresponding average 
velocity of the particles outside of the driven channel, 
(Vf,), versus Fd for a system with a monodisperse trian- 
gular crystalline arrangement of colloids. The different 
curves indicate the effect of changing the interaction co- 
efficient q. For low Fd, all the curves for both (Vc) and 
(Vfc) increase linearly and (Vc) — (Vb), corresponding to 
the elastic flow regime where the particles in the driven 
region are locked with the bulk particles. Here there is no 
tearing or rearrangements, so all the particles keep their 
same neighbors and move together. At higher drives, a 
transition occurs to a regime where (Vb) decreases and 
(Vc) increases with increasing Fd. This occurs when the 
particles in the driven channel partially decouple from 
the bulk particles and are able to move past them. The 
bulk particles continue to move due to the coupling to the 
driven particles; however, as Fd increases, this coupling 
becomes weaker and (Vb) decreases. The inset in Fig.[2lja) 
shows the force Fc at which the decoupling occurs as a 
function of q. As the strength of the particle-particle in- 
teractions increases, a larger force must be applied before 
decoupling can occur. 

The velocity-force curves in Fig. [5] are very similar to 
those found in other systems that exhibit decoupling and 
drag effects. One of the best known examples is the Gi- 



0.0015 

A 

>°o.ooi 

V 

0.0005 
0, 



(a) ' 
0.04* 

*0.02 


L ' l' ' 1 1 
5 2 ^2.5 




.(b) 



'o 0.01 0.02 0.03 0.04 0.05 



Fd 

FIG. 3: A monodisperse system of colloids with q — 1.0 for 
varied colloid density, (a) (Vc) vs Fd- (b) (Vb) vs Fd- The 
colloid lattice constant a = 3.0, 2.5, 2.3, 2.1, 2.0, 1.9, 1.7, 
1.6, and 1.5, from left to right. Samples with larger a are less 
dense and decouple at a lower driving force. The inset in (a) 
shows i^c vs a. 



aever transformer geometry for coupled superconducting 
layers in the presence of a magnetic field^^i^. Each quan- 
tized magnetic field line must pass through both layers, 
providing a magnetic coupling between the layers. When 
only one layer is driven by a current, the voltage response, 
which is proportional to the superconducting vortex ve- 
locity, is identical in both layers, indicating that the vor- 
tices in the two layers are locked together. As the drive 
increases, there is a decoupling transition between the 
layers associated with a decrease in the effective drag on 
the driven layer. This causes the velocity of the vortices 
in the driven layer to increase while the velocity of the 
vortices in the secondary or undriven layer drops. In this 
superconducting system, the interaction between vortices 
in neighboring layers is attractive, while in the colloidal 
system, the interactions between colloids in the driven 
and undriven regions are repulsive. The response of re- 
pulsively interacting particles in coupled one-dimensional 
wires has been studied for classical electrons^ and par- 
ticles with Yukawa interactions^!^. In these systems, a 
commensurate state can form when the number of par- 
ticles in each layer is the same, producing a well defined 
coupling-decoupling transition when one layer is driven. 
In the case of our driven line of particles, the system 
can be viewed as a single driven one-dimensional layer 
interacting with an array of non-driven layers. At the 
decoupling transition, the particles outside of the driven 
region remain locked together so that a very small shear 
band occurs only at the driven line. For more disordered 
systems, a much broader shear band forms. 

In Fig. [3] we show (Vc) and (Vb) vs Fd for a monodis- 
perse sample with fixed q — 1.0 but varied colloid den- 
sity, measured in terms of the lattice constant a. The 
inset of Fig. E^a) shows the dependence of the decou- 



FIG. 4: A bidisperse colloidal system where half the particles 
have charge qi and the other half have charge q2 = 1.0. (a,b,c) 
Dots: Particle positions. Dark (blue) particles have charge qi 
and light (red) particles have charge q2. Lines indicate the 
trajectories followed by the particles over a fixed time interval. 
(d,e,f) Voronoi construction. White polygons correspond to 
particles with six neighbors, dark blue polygons to particles 
with five neighbors, and light red polygons to particles with 
seven neighbors. (a,d) At qi/q2 = 0.7 and Fd — 0.05, the 
system is mostly ordered. There are some dislocations near 
the driven region due to the fact that the driven particles 
are moving faster than the bulk particles. These dislocations 
all have their Burgers vectors aligned in the same direction. 
(b,e) At qi/q2 = 1.6, the system is partially disordered in the 
bulk and there is a mixing of the trajectories as the shear 
band widens. (c,f) At qi/q2 ~ 2.4, the system is strongly 
disordered and the width of the shear band is increased. 

pling force Fc on a. We observe a trend similar to that 
of changing q, where the samples with lower density have 
weaker particle-particle interactions and a lower decou- 
pling threshold. The decoupling transition should also 
depend on the commensurability ratio of the particles in 
the driven channel with the surrounding media. If the 
density in the driven channel is higher or lower than that 
of the bulk, the decoupling force will be depressed rela- 
tive to the case we consider here where the two densities 
are equal. This is a result of the localized incommensura- 
tions that separate the two regions of different densities. 
The incommensurations depin below the bulk decoupling 
transition and shift the transition to lower F^. Studies 
of one-dimensional coupled channels also found that the 
decoupling transition drops to lower drives at incommen- 
surate channel filling ratios^. 



IV. BIDISPERSE SYSTEMS 

We next consider a bidisperse system where half the 
particles have an interaction coefficient of qi and the 
other half have an interaction coefficient of g2- In 
Fig. |ll^a,b,c) we plot the particle positions and trajec- 
tories during a fixed period of time, and in Fig. |31[d,e,f) 



FIG. 5: (Vc) vs Fd for bidisperse systems with §2 = 1.0. (a) 
qi/q2 = 0.1, 0.2, 0.3, 0.5, 0.7, 0.8, 0.9, 1.0, and 1.1, from left to 
right. The system is disordered for qi/q2 < 0.6. (b) (?i/g2 ~ 
1.3 (center right), 1.4 (upper right), and 1.8 (lower right). Fc 
is highest for gi/(/2 ~ 1.3 and lowest for qi/q2 = 1-8. The 
system is disordered for qi/q2 = 1.8, producing the enhanced 
fluctuations in the velocities. There is also a crossing of the 
velocity-force curves, where (Vc) at qi/q2 = 1.8 is lower in 
the moving phase than (Vc) of the ordered state due to the 
widened shear band in the moving state in the disordered 
system, (c) The disordered states at qi/q2 ~ 2.2, 2.4, 2.6, 2.8 
and 3.0, from top to bottom. Here (Vc) for fixed Fd decreases 
with increasing qi/q2- 



we show the Voronoi tessellations of the particle posi- 
tions at one instant. For qi/q2 = 0.7 and Fd — 0.05, in 
Fig. |l{a,d), the system is crystalline with sixfold order- 
ing in the bulk. The topological defects are concentrated 
near the driven line since the system is in the decou- 
pled phase where the particles in the driven line move 
faster than the bulk particles. There is also no trans- 
verse diffusion in the system, as indicated by the nearly 
one-dimensional, non-mixing trajectories of the particles. 
In Fig. |lUb,e) at (71/(72 = 1-6, the bulk is disordered 
with a proliferation of 5-7 paired defects. The trajec- 
tories show that there is strong mixing near the driven 
line, with the most strongly disordered region concen- 
trated near the driven line. Fig. |lljc,f) shows that for 
(71/(72 = 2.4, the system is even more disordered, and 
contains an increased number of 5-7 defects. In addition, 
the region of strong mixing, denoted by the region with 
crossing trajectories, is now wider in the y-direction. For 
0.6 < qi/q2 < 1.6, the bulk is ordered. For (71/(72 < 0.6 
the system disorders and the trajectories are similar to 
those shown in Fig. |H[b) . 

The different phases can also be identified via changes 
in the features of the velocity-force curves. In Fig. ^a.) 
we plot (Vc) vs F]j for samples with (71/(72 = 0.1 to 
1.1. As (71/(72 increases, the decouphng force Fc in- 
creases. The concavity of the velocity-force curves can 
be fit to (Vc) (X {F - FcY, where /3 > 1.0 for the dis- 
ordered systems with (71/(72 < 0.6 and (3 < 1.0 in the 
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FIG. 6: (a) (Vb) vs Fd — Fc for a bidisperse system with 
q2 ~ 1.0, where Fc is the force at decoupling or depinning. (a) 
qi/q2 = 2.6, 2.8, and 3.0, from top right to bottom right. The 
dashed hnes are power law fits (vertically offset for clarity) 
with exponents of /3 = 1.26, 1.29, and 1.32, respectively, (b) 
The same for the ordered regime of qi = 0.7, 0.9, and 1.0, 
from bottom right to top right, with dashed lines indicating 
power law fits (vertically offset for clarity) of /3 = 0.746, 0.84, 
and 0.8, respectively. 



ordered systems. This is shown more clearly in Fig. [5{b) 
where we plot vs Fd for systems with (71/(72 — 1.3, 
1.4, and 1.8. For the disordered case of (71/(72 — 1.8, 
there are stronger fluctuations in (14). Fc is highest for 
Qi/Q2 = 1-3 but is lower in the disordered regime for 
Q1/Q2 — 1-8. A crossing of the velocity-force curves oc- 
curs since Fc for the (71/(72 = 1-8 system is lower than 
for the (71/(72 = 1.3 system, but at higher drives (Vc) 
is lower in the (71/(72 = 1-8 sample. This occurs since 
the ordered state at (71/(72 — 1.3 produces a sharp shear 
band, while the disordered state at 91/(72 = 1.8 has a wide 
shear band region similar to that illustrated in Fig.|4l[b), 
meaning that more particles are being dragged by the 
line particles, causing the driven particles to move more 
slowly. 

Figure^c) shows {Vc) versus Fd for 2.2 < (71/^2 < 3.0, 
where the system remains disordered with large fluctua- 
tions in (Vc). For a fixed drive, (Vc) drops with increas- 
ing qi/q2- Here the velocity-force curves can be fit to the 
form {Vc) - Vdc K {Fd - Fc)^ with /3 « 1.3, where Vdc is 
the velocity at Fc. This fit is illustrated in Fig. IHl^a) for 
(71/92 = 2.6, 2.8, and 3.0. A power law fit of this type 
is similar to that observed for a single particle moving 
through a disordered medium when the driven particle 
creates a large amount of distortion, including particle 
rearrangements^!^. Similar scaling is found at depin- 
ning for assemblies of particles that undergo plastic flow 
upon depinning; such tearing behavior generates an ex- 
ponent ranging from (3 = 1.25 to /3 = 2.022iii. Near an 
elastic depinning transition, the velocity-force curve can 
also be fit with a power law but with a smaller exponent 
/3 < 1. 9 44,45 ^ Pqj. samples in the ordered regime, we can 
fit (Vc) — Vdc to a similar scaling form with j3 « 0.8, as 
shown in Fig.[ni[b) for (71/(72 = 0.7, 0.9, and 1.0. This indi- 
cates that the decoupling transition in the ordered phase 




FIG. 7: (a,b) Velocity profiles Vx vs y for a bidisperse system 
with g2 = 1.0. The driven region is centered at yd = 23. (a) 
In a sample with 91/92 = 1.0 at Fd = 1.0, a sharp shear 
band forms, (b) Samples with 91/92 = 1.8, 2.0, 2.2, 2.4, 
2.6, 2.8, and 3.0, from top center to bottom center. Vx in 
the driven region drops with increasing 91/92, while Vx in 
the bulk increases with increasing 91/92 as the shear band 
widens, (c) Vx in only the top half of the sample vs y ~ yd 
for samples with qi/q2 = 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, and 3.0, 
from top to bottom. Dotted line is a power law fit for the 
91/92 ~ 2.2 system with an exponent of —1.4 in the region 
of strong mixing. At large y — yd, the bulk becomes locked 
again and moves as a solid. 



resembles the elastic depinning of a one-dimensional cou- 
pled chain of particles from a periodic substrate created 
by the periodic ordering of the bulk particles. We find 
similar scalings for the other fillings we have considered, 
where for the ordered systems /? < 1.0 and for the disor- 
dered systems /? > 1.0. 

The shear banding effect can be better seen by ex- 
amining profiles of the average x velocity Vx taken at 
different points along the y-direction. Here, Vx{y) = 
E (Rz • X) [e (i?^y - 27 - (5y /2) e (y + ^ y /2 - i?^ ] where Rl 
is the y coordinate of particle i and 6y is the width of the 
averaging region. Figure [Tja) shows the velocity profiles 
Vx{y) for (71/(72 = 1.0 at Fd = 1.0. A sharp spike in Vx 
appears at y = yd, the center of the driven region, while 
Vx ^ in the bulk undriven portion of the sample. This 
type of profile is observed for all samples where the sys- 
tem remains ordered in the decoupled phase. In Fig. [Tfb) 
we plot Vx{y) for disordered samples with (71/(72 = 1.8, 
2.0, 2.2, 2.4, 2.6, 2.8, and 3.0. In all cases, Vx is maxi- 
mum aX y = yd and falls off for larger \y — yd\- As (71/92 
increases, Vx at y — yd decreases while in the bulk Vx 
increases. This coincides with the decrease in {Vc) for 
increasing 91 / 92 shown in Fig. EJc) . The increase of the 
response in the bulk occurs when a larger number of par- 
ticles are dragged by the driven line due to the increased 
particle-particle interaction strength. In Fig.[7l^c) we plot 
Vx for only the upper half of the sample versus y — yd on 
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FIG. 8: A bidisperse sample with q2 = 1.0. (a) Fc vs qxjqi. 
A peak appears in the ordered region centered close to qxjqi = 
1.0. (b) (Vb) vs gi/(/2 at Fd = 0.05. The velocity drops in the 
disordered regime qxjqi > 1.6. The dashed line indicates a 
fit to 1/(51/92) in the disordered regime, (c) The standard 
deviation of the fluctuations in the bulk velocity 5Vb vs qi jqi 
for Fi = 0.05. (d) SVt vs qi/q2 for Fa = 0.1. 

a log-log scale, to better shovif vifhere the shear banding is 
occurring. The dashed line in Fig. [TJ^c) is a power law fit 
to the data from the qi/q2 — 2.2 sample in the shear band 
region, and has an exponent of —1.4. From the trajec- 
tory images in Fig. |4jc) , it is clear that there is a region 
of strong mixing near the driven line and a region of less 
mixing at larger values of y — y^, correlated with a satura- 
tion of Vx in Fig.[7l^c). As qi/q2 increases, the saturation 
region shifts out to larger values oi y — i/d- This results 
shows that the system has liquid-like behavior near the 
shear band region, but acts more like a disordered elastic 
solid away from the driven region. 

We summarize the behavior of Fc vs qi / q2 for the bidis- 
perse samples in Fig. [SJ^a). In the disordered regime 
0.1 < (71/(72 < 0.6, Fc increases with increasing (71/(72. 
Above (71/92 = 0.6, when the system enters the ordered 
regime, F^ increases more rapidly and the crystalline or- 
der that appears near (71/(72 = 1-0 enhances the cou- 
pling transition, as indicated by the peak in Fc near 
(71/(72 = 1-0. As the system becomes disordered away 
from (71/(72 = 1.0, weak or defected spots appear that 
lower the decoupling transition and reduce Fc. A min- 
imum in Fc appears near (71/(72 = 1.8 in the disordered 
region, while Fc increases again for larger 51/(72 due to 
the increasing strength of the particle-particle interac- 
tions, which also cause the system to depin plastically. 
The inset of Fig. ^a) shows that Fc in the monodisperse 
system also increases with increasing q. 

In Fig. IHlJb) we plot (Vf,) vs (71/92 in the bidisperse 
samples at Fd = 0.05. Here the bulk velocity decreases 
at a moderate rate with increasing (71/(72 until the sys- 
tem becomes disordered for (71/(72 > 1.6 and large shear 
bands form, leading to a more rapid drop in (Vt) with 




FIG. 9: (a,b) Power spectra S{i^) for bidisperse samples with 
q2 = 1.0 at Fd — 0.05. (a) In ordered systems, (71/(72 = 0.6 
(black), 0.8 (red), 1.0 (green), and 1.4 (blue), there is a narrow 
band noise signature, (b) In disordered systems, gi/g2 = 2.4 
(black), 2.6 (red), 2.8 (green), and 3.0 (blue), a broad band 
noise signal occurs. The dashed line is a power law fit with 
an exponent of —1.5. (c) The noise power So vs q\/q2- At 
the transition into the disordered regime. So jumps to a much 
higher value. 



increasing qxjqi- The dashed line is a fit to l/{qi/q2) in 
the disordered regime. Just before the system disorders, 
for 1.4 < (71/92 < 1-6, there is a slight increase in (Vf,) 
that is correlated with a drop in the standard deviation of 
the velocity fluctuations (514, as shown in Fig. El^c). For 
Fd ~ 0.05, Fig. Wis) indicates that 5Vb reaches a maxi- 
mum for 91/92 = 1-1, the same value at which there is a 
peak in Fc in Fig. HJa). At this point, when the veloc- 
ity oscillations are the largest, the particles are moving 
in a one-dimensional ordered pattern and each particle 
exhibits an oscillating velocity component due to the ef- 
fective periodic potential created by the ordered particles 
in the bulk. There is a dip in 5Vb at the transition from 
the ordered to the disordered region at 91/92 = 1-6; for 
91/92 > 1-6, 5Vb decreases with increasing 91/92. At the 
disordering transition, the particles no longer move to- 
gether but have different sliding velocities, resulting in 
a cancellation of the large oscillations that appeared in 
the ordered phase. For larger 91/92, deeper into the dis- 
ordered regime, the velocity fluctuations decrease with 
increasing 91 /92 due to a decrease in the amount of mo- 
tion in the system as the flxed value of Fd gets closer to 
Fc^ which increases with increasing 91/92. At Fd — 0.1, 
Fig. [DJd) shows that the same general trends in 5Vb per- 
sist; however, in this case the velocity fluctuations for 
larger 91/92 increase with increasing 91 since the higher 
drive causes a large shear banding effect. There is a sat- 
uration of 5Vb near 91/92 = 3.0, and we expect that the 
velocity fluctuations will decrease for higher 91/92 as Fc 
increases and approaches the flxed value Fd = 0.1. These 
results indicate that the order-disorder transition can be 
detected by measuring the velocity fluctuations. 

The velocity fluctuations can also be characterized us- 
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ing the power spectrum S{iy) = \ J Vb(t)e~'^'^^'^^dt\'^ . In 
the ordered regime, the noise fluctuations have a narrow 
band noise feature indicating that there is a characteris- 
tic frequency. Similar narrow band noise occurs for par- 
ticles moving over a periodic substrate, where the char- 
acteristic frequency is determined by the rate at which 
the particle travels from one substrate minimum to the 
next^^. In Fig. HJa) we plot S{i') in the ordered regime 
for samples with qi/q2 = 0.6, 0.8, 1.0, and 1.4 for a fixed 
drive of = 0.05. The peaks in the spectra indicate the 
presence of narrow band noise. In the disordered regime, 
we find broad band noise or a 1//" noise signature, as 
shown in Fig.HJb) for qi/q2 = 2.4, 2.6, 2.8, and 3.0. Here 
a = —1.5. Similar broad band noise has been observed 
for dragging a single particle through granular media at 
the jamming transitions and for plastic depinning of par- 
ticles on disordered substrates^. We can also analyze the 
noise power So = J^^ Slv) for a fixed octave, as shown 
in Fig. 9(c) for i^i = 10 and 1^2 = 100. Here Sq is small 
in the ordered regime 0.6 < qi/q2 < 1-6 and undergoes a 
pronounced increase to a maximum near the onset of the 
disordered regime. For qi/q2 < 0.6, the system is also 
disordered and S{i') shows broad band noise; however, 
the noise power remains very low for these small values 
of qi/q2- 



A. Dynamic phases 

For bidisperse samples with large qi/q2, we find an ad- 
ditional dynamical phase at high drives where the shear 
band region becomes strongly localized again and the 
bulk particles lock together, similar to the behavior found 
in the ordered regime. In Fig. [TUf a) we plot (Vc) versus 
Fd for a bidisperse sample with qi/q2 = 1.6. For low Fd, 
the system is crystalline and strongly coupled, so that 
all the particles move together. As Fd increases, the sys- 
tem becomes disordered and produces a large shear band 
(illustrated in Fig. [TUTc) for Fd — 0.05) associated with 
large velocity fluctuations. At higher drives, however, 
there is a transition to a strongly shear localized state 
where only a single line of particles are moving while the 
bulk particles lock together and moving at a very slow 
velocity, as shown in Fig. fTUTd) for Fd = 0.09. The onset 
of this phase is accompanied by a drop in the velocity 
fluctuations in Fig.fTUTa). In Fig.lTUTb) we plot (Vf,) ver- 
sus Fd for the same system, showing that at the shear 
localization transition, the bulk velocity drops to a value 
slightly above zero. At the transition, (Vc) increases since 
the drag from the bulk particles on the particles in the 
driven channel is reduced. The shear localization occurs 
when the particles in the driven line are moving suffi- 
ciently rapidly that they can no longer couple effectively 
to the bulk particles. This transition is a general feature 
that occurs in the disordered state. Since the transition 
is dominated by fiuctuations, the drive at which it oc- 
curs can vary significantly from one sample realization 
to another; however, on average the transition occurs at 





FIG. 10: (a) (Vc) vs Fd for a bidisperse system with qi/q2 = 
1.6 and 52 ~ 1.0. At low drives, the system is crystalline 
and moves as a solid unit. At decoupling, the system enters 
a disordered regime, while at higher drives Fd > 0.075, there 
is a transition to a state with strongly localized shear and 
reduced velocity fluctuations, (b) (Vi,) vs Fd for the same 
system. The bulk velocity drops almost to zero at the onset 
of the shear localization transition. (c,d) Particle positions 
(filled circles) and trajectories (lines) in the same system for 
(c) the disordered decoupled state at Fd = 0.05 and (d) the 
strongly localized shear state at Fd — 0.09. Arrows in (a) and 
(b) indicate the drives at which the images in (c) and (d) were 
obtained. 

higher values of Fd for increasing qi/q2 or decreasing a. 
The shear localized state has small velocity fiuctuations 
with low noise that is white (a = 0). 

V. QUENCHED DISORDER AND SHEAR 
LOCALIZATION 

We next examine the effects of adding quenched dis- 
order or pinning sites in the bulk region. In Fig. 11 we 
plot the particle trajectories in the disordered regime for 
a bidisperse sample with qi/q2 — 1.8 at Fd ~ 0.05. The 
pinning sites are placed an average distance dp from the 
driven region. Each pin captures a single colloid and the 
pinning force is sufficiently strong that the particles do 
not depin over the range of Fd we consider. The repul- 
sive nature of the particle-particle interactions prevents 
unpinned particles from closely approaching pinned par- 
ticles, leading to a reduction in the density of particle 
trajectories in the vicinity of each pinned particle. For a 
small number of pinning sites Np = 4, shown in Fig. lllf a) 
for a sample with dp = 3a, the overall particle trajectories 
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FIG. 12: (a) (Vc) vs Fd for a bidisperse system at gi/(/2 = 1-8 
and 52 ~ 1-0 with quenched disorder where A'j, = 16 and 
dp = 2a. The pinned, coupled disordered, and high drive de- 
coupled regimes are labeled, (b) Dynamic phase diagram for 
Fd vs A'j, in bidisperse samples with qi/q2 = 1-8. The pinned, 
coupled disordered, and high drive decoupled regimes are la- 
beled. As Np increases, the transition between the coupled 
and decoupled phases drops to lower values of Fd- 



FIG. 11: The particle positions (filled circles) and trajecto- 
ries (lines) at a constant drive of Fd = 0.05 for a bidisperse 
system at qi/q2 = 1-8 and q2 = 1.0 with quenched disorder 
added in the form of localized pinning sites placed at an av- 
erage distance of dp from the driven line. The open circles 
are the locations of the pinning sites that each capture one 
colloid, (a) A'j, — 4 and dp = 3a. (b) A'p = 20 and dp — 3a. 
As the number of pinning sites increases, the flow becomes 
more localized, (c) Np = 10 and dp = 2a. (d) Np — 10 and 
dp — 18a. 



do not differ significantly from the pin-free system; how- 
ever, as the number of pinning sites increases, the mixing 
region surrounding the driven line is reduced in width, as 
illustrated in Fig. [TlT b) for a sample with Np = 20 and 
dp = 3a. Here, there is no longer any net motion of 
the bulk particles, indicating a complete screening of the 
shear banding effect by the pinning. This result shows 
that quenched disorder can produce strong shear local- 
ization. 

In Fig. [T^ a) we plot (Vc) versus Fd for a bidisperse 
pinned system with qi/q2 — 1-8, Np — 16, and dp — 2a 
where we observe a pinned regime, a coupled disordered 
flow regime, and a decoupled regime for high drive simi- 
lar to that found for the pin-free system in Fig.[TUl In the 
absence of pinning, at low drives the particles all move 
as a locked solid; however, when pinning is added to the 
system, there is a true pinned phase at low drive where 
particle motion does not occur. In Fig. [T^ b) we plot 
the dynamic phase diagram for Fd versus the number of 
pinning sites A'j,. The transitions between the different 
phases show strong fluctuations from one realization to 
another, and in Fig. [T^b) we show the transition line 
averaged over several different disorder realizations. The 
transition between the coupled and decoupled regimes 



drops to lower Fd with increasing Np. For a fixed drive, 
we find the interesting effect that the onset of the decou- 
pled regime is correlated with an increase in the velocity 
of the particles in the driven line. Thus, by increasing the 
number of pinning sites, it is possible to cause the driven 
particles to move at a higher velocity. This effect arises 
due to the effective screening of the driven particles from 
the bulk particles by the pinned particles. A similar ef- 
fect was observed for driving a single particle through a 
background of other particles and pinning sites, where for 
high drives, a decoupling transition between the driven 
and background particles leads to an increase in the ve- 
locity of the driven particle. In Fig. [TIT c.d) we show the 
particle positions and trajectories for samples where the 
number of pinning sites is fixed at Np = 10 for differ- 
ent pin spacings of dp — 2a [Fig. Illf c)] and dp — 18a 
[Fig. rnT d)]. Here the motion is more localized for the 
smaller value of dp . 

In Fig.fTHTa) we show Vx{y) for bidisperse samples with 
91/92 = 1-8, dp — 2a, and Fd — 0.024 for different num- 
bers of pinning sites ranging from A'j, = to Np = 20. 
For the pin-free system with Np — 0, Vx far from the 
driven line reaches a finite value since all of the particles 
in the system move. As the number of pinning sites in- 
creases, Vx decreases at all length scales, with Vx in the 
driven line region decreasing more slowly than Vx in the 
region far from the driven line. For a sufficiently large 
number of pinning sites, Vx drops to zero for distances 
of 6a or greater away from the driven line, indicating a 
complete screening of the shear band by the pinned par- 
ticles. In Fig.fTWb) we show Vx for the same system with 
fixed A'j, = 10 and varied dp. For increasing dp, Vx in the 
bulk gradually increases. 
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FIG. 13: Velocity profiles Vx vs y for bidisperse samples 
with qi/q2 ~ 1-8 and q2 = 1.0 at Fd = 0.024. (a) dp = 2a and 
Np = 0, 4, 6, 8, 10, 12, 14, 16, 18, and 20, from top center 
to bottom center. Here Vx drops with increasing Np at all 
length scales, (b) Np — 10 and dp = 2a, 4a, 6a, 8a, 10a, 12a, 
14a, 16a, and 18o, from bottom center to top center. For the 
smallest dp, the shear band is strongly localized. 



VI. THERMAL EFFECTS 

To study thermal effects, we focus on the monodisperse 
system with qi/q2 = 1-0 under fixed drive Fd = 0.05, 
which forms a crystalhne state. In Fig. [T3] we plot the 
particle positions and trajectories as well as Voronoi con- 
structions of snapshots of the particle positions for in- 
creasing T. For low temperatures T < 0.6, the system 
remains ordered and the particle trajectories are one- 
dimensional. At T = 0.6, Fig. [HTa.d) shows that the 
trajectories near the driven line begin to mix, produc- 
ing 5-7 paired defects, while the bulk particles remain 
ordered. At T = 1.2 in Fig. \MUo,e), the number of 5-7 
defect pairs has increased and some pairs begin to mi- 
grate from the driven line into the bulk region, although 
the bulk remains mostly ordered. For T = 1.56 the bulk 
begins to disorder and a widened shear band appears. 
Figure [Tif cf) shows that at T = 1.8 the bulk is strongly 
disordered. 

In Fig. [TSla'). the plot of (14) versus T shows that the 
velocity drops when the dislocations begin to migrate out 
from the driven line region. There is a minimum in (Vc) 
at the bulk melting temperature. Figure ITST b) shows the 
corresponding fraction of six-fold coordinated particles, 
Pg = X^iii ^i^i where Zi is the coordination number 
of particle i. For a perfect triangular lattice, Pq = 1.0. 
We find that the drop in Pq, indicating bulk melting, oc- 
curs at a higher temperature than the drop in (Vc). In 
previous numerical studies of single driven probe parti- 
cles moving through a crystalline background, it was also 
found that the velocity drops when local dislocations near 
the driven probe particle occur, which typically begins 
for temperatures below the bulk melting temperaturei. 



y 




(a) X (b) X (c) X 




FIG. 14: A monodisperse pin-free system with qi/q2 ~ 1.0, 
q2 — 1.0, and Fd = 0.05 at different temperatures. (a,b,c) 
Dots: Particle positions. Lines: the trajectories followed by 
the particles over a fixed time interval. (d,e,f) Voronoi con- 
struction. White polygons correspond to particles with six 
neighbors, dark blue polygons to particles with five neigh- 
bors, and light red polygons to particles with seven neighbors. 
(a,d) T — 0.6. The dislocations are primarily located along 
the driven line with their Burgers vectors aligned in the same 
direction. (b,e) T — 1.2. The dislocations have begun to mi- 
grate into the bulk. (c,f) T = 1.8, above the bulk melting 
temperature. The fiow is disordered and topological defects 
have proliferated. 




T 



FIG. 15: (a) (Vc) vs T for a pin-free monodisperse system 
with gi/<j2 = 1.0 and 172 = 1.0 at Fd = 0.05. (b) The fraction 
of six-fold coordinated particles P(, vs T. The system melts 
near T = 1.5 and the velocity of the particles in the driven line 
drops at temperatures just below the bulk melting transition. 

At much higher temperatures, the velocities gradually 
increase again. This result shows that the driven line 
geometry can also be used to examine thermal melting 
properties. 



VII. SUMMARY 

We study the effects of driving a quasi-one-dimensional 
region of particles through a two-dimensional system of 
Yukawa interacting particles. In a monodisperse system 
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a crystalline state forms, and as a function of increas- 
ing external drive there is a well defined transition from 
elastic flow, where all the particles move together, to a 
decoupled flow, where the particles in the driven region 
decouple from the bulk particles and the bulk particles 
remain in a locked crystalline state. The properties of 
this elastic to decoupled transition are similar to those 
found in studies of layered systems where only one layer 
is driven and the other layer is dragged, such as a trans- 
former geometry for bilayer superconductors or bilayer 
Wigner crystals. For bidisperse systems, the bulk be- 
comes disordered, and the driven line produces a local 
shear band and velocity gradient. The bulk particles near 
the driven line can be dragged with the driven particles, 
and the average velocity of the particles in the driven line 
decreases as the width of the shear band region increases. 
In the disordered regime, for increasing drive we identify 
another decoupling transition where the bulk particles 
become locked, producing an elastic disordered solid, co- 
inciding with the shear band becoming very sharp. We 
show how the decoupling force, noise fluctuations, and 
average velocities can be correlated with bulk disorder- 
ing transitions. We also consider the effects of adding 



quenched disorder or pinning, and find that the shear 
band region becomes increasingly localized for increas- 
ing pinning density. This results in an interesting effect 
where the velocity of the driven particles can be increased 
by adding more pinning sites to the system. Such sys- 
tems could be realized experimentally using colloids or 
dusty plasmas with an optical drive applied along a one- 
dimensional channel, such as a laser applied along the 
edge of the sample. Variations of this geometry could 
also be constructed in superconducting vortex systems 
by embedding a single weak pinning channel in a bulk 
sample. 
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